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Abstract. Stability of radial and nonradial oscillations of massive supergiants is discussed. 
The kappa-mechanism and strange-mode instability excite oscillations having various periods in 
wide ranges of the upper part of the HR diagram. In addition, in very luminous (log L/Lq > 
5.9) models, monotonously unstable modes exist, which probably indicates the occurrence of 
optically thick winds. The instability boundary is not far from the Humphreys-Davidson limit. 
Furthermore, it is found that there exist low-degree (£ = 1,2) oscillatory convection modes 
associated with the Fe-opacity peak convection zone, and they can emerge to the stellar surface 
so that they are very likely observable in a considerable range in the HR diagram. The convection 
modes have periods similar to g-modes, and their growth-times are comparable to the periods. 
Theoretical predictions are compared with some of the supergiant variables. 
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1. Introduction 

Light and velocity variations on various time-scales are common in very luminous 
stars (e.g., van Genderen 1989, van Leeuwen et al 1998). Those variations are caused 
by various kinds of instabilities. Here we discuss mainly the cause of microvariations of 
massive supergiants based on linear stability analyses applied to evolutionary models 
of massive stars. The evolution models were calculated by a Henyey-type code using 
OPAL opacity tables (Iglesias & Rogers 1996). Wind mass-loss is included for the evo- 
lutionary models for Mi > 30M Q (Mi=initial mass) based on the mass- loss rates of 



Vink et al. (2001) Linear stability analyses were performed using the methods given in 



Saio, Winget & Robinson (1983) for radial modes and Saio & Cox (1980) for nonradial 
modes, where the outer mechanical boundary condition was modified to <5-P gas — > (<5P gas 
means the Lagrangean perturbation of gas pressure) taking into account the fact that 
radiation pressure is dominant near the outer boundary. It should be noted, however, 
that the effect of winds on radial and nonradial oscillations is not included because the 
effect is not well understood yet. 



2. Stability of radial modes 

Figure [1] shows instability boundaries of spherical symmetric modes and selected evolu- 
tionary tracks. Short-dashed line indicates the instability boundary for low-order modes. 
The nearly vertical "finger" (4.4 > logT c ff > 4.3) is the well-known j3 Cep instabil- 
ity strip, in which low-order radial modes and nonradial p-modes are excited by the 
kappa-mechanism at the Fe-peak of opacity around T ~ 2 x 10 5 K. At the luminous 
part of the instability strip, the cool-side boundary bents to become horizontal around 
log(L/L ) ~ 5.8. This is due to the strange-mode instability which occurs in models 
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Figure 1. Instability boundaries of radial modes and selected evolutionary tracks. Short- and 
long-dashed lines indicate stability boundaries for low-order and relatively high-order radial 
modes, respectively. Dotted line indicates the boundary for monotonously unstable radial modes. 
Numbers along evolutionary tracks indicate initial masses in solar units. The instability bound- 
aries for radial modes also represent approximately the boundaries for low-degree (£ < 2) non- 
radial p-modes. 

with sufficiently high luminosity to mass ratios as L/M > 10 4 L q /Mq. The proper- 



ties of the strange modes have been investigated by e.g., Glatzel & Kiriakidis (1993) 



Glatzel (1994)1 |Saio, Baker fc Gautschy (1998)| (see also Saio 2008). 

Nearly vertical part of the boundary around log T c g ~ 3.8 indicates the well-known blue 
edge of the Cepheid instability strip. (Red edge is not obtained because the perturbation 
of convective flux is neglected in this analysis.) 

In the vertical narrow region indicated by long-dashed line around logT ff ~ 3.9 — 3.8, 
relatively high-order radial modes are excited around hydrogen ionization zone. (Low- 
degree nonradial modes with similar frequencies are also excited.) The amplitude of these 
modes are extremely confined to the outermost layers above the hydrogen ionization zone. 
Since these modes exist even under the NAR approximation where thermal-time is set 
to be zero (Gautschy & Glatzel 1990), they may be classified as strange modes. These 
modes have got little attention so far (cf. Gautschy 2009). 

Dotted line in Fig. Q] shows the boundary above which monotonously unstable modes 
exist. The growth times of these modes are much shorter than the timescale of evolution. 
The presence of such monotonously unstable modes have not been recognized before. Such 
a mode probably corresponds to the presence of an optically thick wind as investigated 



by Kato & Iben (1992)1 for WR stars. It is interesting to note that the boundary is not 



far from the Humphreys-Davidson limit (Humphreys & Davidson 1979). 



3. Stability of nonradial modes 

The instability ranges of low-degree high-order g-modes are shown in Fig. [2] by solid 
and dotted lines for I = 2 and I = 1 modes, respectively. These modes are excited 
by the kappa-mechanism at the Fe opacity peak. They are responsible for the long- 
period variations in slowly pulsating B (SPB) stars. Such g-modes can be excited even in 
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Figure 2. Ranges where nonradial oscillations are expected to be detected. The instability 
boundaries for SPB-type low-degree high-order g-modes are shown by solid (t = 2) and dotted 
(t = 1) lines. The ranges where oscillatory convection modes should be observable are indicated 
by long- and short-dashed lines for 1=1 and i = 2, respectively. 



supergiants (SPBsg) because a shell convection zone associated with hydrogen burning 
reflects some g-modes and hence suppresses dissipation otherwise expected in the core 
(Saio et al 2006). 

The long dashed and short dashed lines in Fig. [2] indicate ranges where oscillatory 
convection modes of i = 1 and £ — 2, respectively, are expected to be observable. It 
is well known that linear convection modes are monotonously (dynamically) unstable 
in adiabatic analyses. Shibahashi & Osaki (1981) found, however, that the high-degree 
(£ ^ 10) convection modes become overstable (oscillatory) when the nonadiabatic effect 
is included in luminous (L/Lq = 10 5 ) models hotter than the cepheid instability strip. 

In our massive star models, it is found that low-degree (£ ^ 2) oscillatory convection 
modes exist associated with the Fe-opacity peak convection zones, and some of these 
modes are expected to be observable. Since the growth time of a convection mode is 
short (comparable to the period), the mode is expected to have a large amplitude in 
the convective zone. Therefore, the visibility of the oscillatory convective modes can be 
measured by the ratio of the photospheric amplitude to the maximum amplitude in the 
interior (mostly in the convection zone). Assuming that an oscillatory convection mode 
is observable when the ratio is larger than 0.2, the boundaries of the visible ranges are 
shown in Fig. [2] by dashed lines. Oscillatory convection modes are visible in sufficiently 
luminous (log L/L@ > 4.6) B-type stars. Although there are many oscillatory convection 
modes in a star, only one or two modes for a given I are visible because the other modes 
are well confined to the convection zone. Periods of these modes are comparable to g- 
modes much longer than those of radial modes. These oscillatory convection modes might 
be responsible for long-period variations in supergiant stars. 
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Figure 3. Ranges where (quasi-) periodic variabl e stars are expected o n the HR diagram are 
compare d with supergiant var iables analyzed by Lefever et al. (2007) (big filled circles) and 
LBVs by Lamers et al. (1998) (filled squares). Positions at different LBV (or S Dor) phases for 
each star are connected with thin lines; 1 = R 712, 2 = HR Car, 3 = 164 G Sco, 4 = S Dor, 5 
= R 1273, 6 = AG Car. 



4. Comparison with supergiant variables 

Observed positions of variable supergiants on the HR diagram are compared with the 
instability and visibility boundaries in Fig. [3] Their periods-T c ff relations are compared 
with theoretical ones in Figs. 0] and [SJ Fig. [3] indicates that all the hotter (logT c ff > 4) 
and luminous (log L/Lq > 4.5) stars are expected to show (quasi-)periodic variations, 
which is consistent with the observational fact that no or at most a very little number of 



stable supergiants exist in a spectral range of 09 - AO as found by van Genderen (1989) 



Big dots in Fig.[3]are relatively less luminous supergiant variables analyzed by Lefever et al. (2007) 



This figure indicates most of them to have masses ranging from ~ 15M Q to ~ 30M Q . 
They are located on the HR diagram in the g-modc instability regions or visible range 
of oscillatory convection modes. Fig. @] compares the periods of these stars as function of 
the effective temperature with theoretical ones of 15M and 3OM models. This figure 
indicates that for most of these stars periods seem consistent to low-degree high-order 
g-modes (SPBsg) or oscillatory convection modes. We note, however, that for the coolest 
three stars the periods are shorter than any of the excited modes. Although these three 
stars are located on the HR diagram in the region where oscillatory convection modes 
should be visible (Fig. [3]), the periods are much shorter than those of the oscillatory 
convection modes. If these effective temperatures are accurate, an unknown excitation 
mechanism might be working in these stars. 

Fig- HI compares theoretical periods of very massive models (Mj = 50M Q and 70M Q ) 
with periods of microvariations of some LVB stars analyzed by Lamers et al. (1998)| each 
of which has different periods and effective temperature depending on the LBV (S Dor) 
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Figure 4. Periods of excited modes and of observable oscillatory convection modes versus 
effective temperature, where open (30Mq) and filled (15M@) circles are radial modes; triangles 
(30Mq) and inverted triangles (15M©) are 1=1 modes; squares (30Mq) and diamonds (15Mq) 
are I — 2 modes. Convection modes are shown by filled symbols for nonradia l modes. Big dots 



show observed periods-T c ff relations of supergiants from Lefever et al. (2007) 



phases. Figs. [5] and [3] indicate that the microvariations of R 712 (1), S Dor (4) and AG Car 
(6) are consistent to the properties of oscillatory convection modes. 

The periods of microvariations of HR Car (2), 164 G Sco (3), and AG Car (5) are 
consistent with periods of strange modes. However, luminosities of HR Car and 164 G Sco 
are too low for the strange modes to exist (Fig. [3]) . Further investigations are needed to 
resolve the discrepancy. 



5. Summary 

We have discussed various instabilities that occur in massive supergiants. Radial modes 
and nonradial p- and g-modes are excited by the kappa-mechanism at the Fe opacity 
bump at T ~ 2 x 10 5 . 

In a star with a very high luminosity to mass ratio of L/M > 10 4 L©/M©, the strange 
mode instability works for radial and nonradial modes. Strange modes seem to be re- 
sponsible for quasi-periodic variations in some of the luminous supergiants. 

In addition, it is found that in very luminous models (logL / X© > 5.9) a monotonously 
unstable radial mode exists, which is probably related to the occurrence of an optically 
thick wind. It is interesting to note that the instability boundary roughly coincides with 
the Humphreys-Davidson limit. 

Furthermore, we found that low-degree (£ = 1,2) oscillatory convection modes exist in 
the convection zones caused by the Fe opacity peak, and that some of them can emerge 
to the stellar surface and hence be observable. The oscillatory convection modes have 
periods of 10 ~ 10 2 days depending of the effective temperature, which are longer than 
those of strange modes. The growth-times are comparable to the periods. They seem to 
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Figure 5. Theoretical periods of excited modes and of observable oscillatory convection modes 
versus effective temperature, where open (70Mq) and filled (50Mq) circles are radial modes; 
triangles (70Mq) and inverted triangles (50Mq) are 1=1 modes; and squares (JQMq) and 
diamonds (50M©) are I = 2 modes. Convection modes are shown by filled symbols for nonradial 
modes. Big filled squares connected with thin lines show observed periods- T c s relations of LBVs 



obtained by Lamers et al. (1998) Periods at different LBV (or S Dor) phases for each star are 
connected with thin lines; 1 = R 712, 2 = HR Car, 3 = 164 G Sco, 4 = S Dor, 5 = R 1273, 6 
= AG Car. 

be consistent with the properties of long-period microvariations in LVB stars (see Saio 
2010 for further discussions). 
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